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ABSTRACT 

Using a sample of 69,919 red giants from the SDSS-III/APOGEE Data Release 12, we measure the 
distribution of stars in the [a/Fe] vs. [Fe/H] plane and the metallicity distribution functions (MDF) 
across an unprecedented volume of the Milky Way disk, with radius 3 < R < 15 kpc and height 
\z\ < 2 kpc. Stars in the inner disk (R < 5 kpc) he along a single track in [a/Fe] vs. [Fe/H] , starting 
with a-enhanced, metal-poor stars and ending at [a/Fe] ^ 0 and [Fe/H]~ +0.4. At larger radii we 
find two distinct sequences in [a/Fe] vs. [Fe/H] space, with a roughly solar-a sequence that spans 
a decade in metallicity and a high-a sequence that merges with the low-a sequence at super-solar 
[Fe/H]. The location of the high-a sequence is nearly constant across the disk, however there are 
very few high-a stars at R > 11 kpc. The peak of the midplane MDF shifts to lower metallicity at 
larger R, reflecting the Galactic metallicity gradient. Most strikingly, the shape of the midplane MDF 
changes systematically with radius, with a negatively skewed distribution at 3 < R < 7 kpc, to a 
roughly Gaussian distribution at the solar annulus, to a positively skewed shape in the outer Galaxy. 

For stars with M > 1 kpc or [a/Fe] > 0.18, the MDF shows little dependence on R. The positive 
skewness of the outer disk MDF may be a signature of radial migration; we show that blurring of 
stellar populations by orbital eccentricities is not enough to explain the reversal of MDF shape but a 
simple model of radial migration can do so. 

Subject headings: Galaxy: abundances, Galaxy :disk, Galaxy: stellar content, Galaxy structure 
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1. INTRODUCTION 

The Milky Way is an excellent testing ground of our 
understanding of Galaxy evolution, due to the ability to 
resolve individual stars and study stellar populations in 
greater detail than in other galaxies. However, under¬ 
standing the composition, structure, and origin of the 
Milky Way disk is still currently one of the outstanding 
questions facing astro nomy, an d there is great debate 
about this topic (e.g., Rix & Bovy 2013). The ability 
to resolve individual stars allows one"~to trace the fossil 
record of the Milky Way across the disk, as the stars con¬ 
tain the chemical footprint of the gas from which they 
formed. 

Observations of stars in the Milky Way have led to the 
discovery of several chemical and kinematic prope rties o f 
the d i sk of the Galaxy, such as the thick disk (e.g., Yoshii 
19821 iGilm ore fe R eid 198 3]), chem i cal abundance gradi- 
Hartkopf &; Yoss 1982; |Cheng et al.| |2012b 

2i 


ents 


e.g. 


Anders et al.lh 


mm. 

2014|), and the GMwar f problem ([van den Berghf. 

Pag e! fc Patchett|1975 ) from the study of the metal 


Hayden et al.||20141 ISchlesinger et al 


[1962 

licity 


distribution funct ion (MPF) of the solar neighborhood 

(jvan den Bergh|[l962l Casagrande et al.||2011 Lee et al 


Oil; Schlesinger et al. 2012p . Much of the previous work 


on the Milky Way disk has focused on the solar neighbor 
hood, or tracer populations (e.g., Cepheid variables, HII 
regions) that span a narrow range in age and number only 
a few hundred objects even in the most thorough studies^ 


The advent of large sc ale surveys such as SEG UE (Yanny 
eTnj|2009l), RAVE J Steinmetz et aLl|2006|), APOGEE 


(Maiewskiet al.|2015|), GA1A-ESO (IGilmore et al.|2012|) , 

U L'O I W I L'O n A I A n * 


and HERMES-GAL AH ( Freeman et al.|2012| ) aim to ex- 


pand observations across the Milky Way and will greatly 
increase the spatial coverage of the Galaxy with large 
numbers of stars. In this paper we use observations from 
the twelfth data release of SDSS-III/APOGEE (Al 


am 


et al. 2015| ) to measure the distribution of stars in the 


a/Fej vs. [Fe/H] plane and the metallicity distribution 
function across the Milky Way galaxy with large numbers 
of stars over the whole radial range of the disk, [a/ H] is 
defined by the elements O, Mg, Si, S, Ca, and Ti changing 
together with solar proportions. Using standard chemi¬ 
cal abundance bracket notation, [a/Fe] is [a/H]-[Fe/H]. 

Different stellar populations can be identified in chem¬ 
ical abundance space, with the a abundance of stars 
separating differing populations. The distribution of 
stars in the [a/Fe] vs. [Fe/H] plane shows two dis- 


Fuhrmann 1998; 

Prochaska et al. 2000 

Rede 

y et al. 2006 

Adibekyan et al. 

2012| Haywood et al. 

2(H3f 

Anders et al. 

2014 

Bensby et al.|2014 Nidever et al.|2014 

Snaith et al. 


across a large range of metallicities, and the other track 
having a high-[a/Fe] ratio at low metallicity that is con¬ 
stant with [Fe/H] until [Fe/H] ~ —0.5, at which point 
there is a knee and the [a/Fe] ratio decreases at a con¬ 
stant rate as a function of [Fe/H] , eventually merging 
with the solar- [a/Fe] track at [Fe/H] ~ 0.2 dex. The 
knee in the high- [a/Fe] sequence is likely caused by the 
delay time for the onset of Type la SNe (SNela): prior 
to formation of the knee, core collapse supernovae (SNII) 
are the primary source of metals in the ISM, while after 
the knee SNela begin to contribute metals, enriching the 


ISM primarily in iron peak elements and lowering the 
[a/Fe] ratio. 

Stars on the [a/Fe] -enhanced track have much la rger 
vertical scale-heights than solar-[a/ Fe] stars (e.g., Lee 


,_,,_,,_,,_ Ml. _ 

et al. 2011 Bovy et al. 2012a|b|c|) and make up the 
stellar populatio ns belonging to the thick disk. Nide- 
ver e t al. ([20141) used the APOGEE Red Clump Gat- 
alog (jBovy et al. 2014|) to analyze the stellar distribu¬ 
tion in the [a/Fe] vs. [Fe/H] plane across the Galac¬ 
tic disk, and found that the high-[a/Fe] (thick disk) se¬ 
quence was similar over the radial range covered in their 
analysis (5 < R < 11 kpc). The thick disk stellar pop¬ 
ulations are in general observed to be more metal-poor, 
a-enhanced, have shorter radial scale-lengths, larger ver- 


in the solar neighborhood (e.g., Bensby et al. 

2003; Al- 

lende Prieto et al. 

2006; Bensby et al. ‘ 

2011i b 

>ovy et al. 

2012c Gheng et al. 

2012a Anders et a 

. 2014), although 

there do exist thic 

k disk stars with sc 

)lar-|a/Fe] abun- 


dances and super-solar metal l icities (|Bensby et al.||2 003 
20051 jAdibekyan et al.||2011 JBensby et al.||2014| |JNicfe^ 


ver et al.||2014; Snaith et al. 2014|). However, the exact 
structure of the disk is still unknown, and it is unclear 
whether the disk is the superposition of multiple com¬ 
ponents (i.e., a thick and thin disk), or if the disk is a 
continuous sequence of stellar pop ulations (e.g., Ivezic 
et al. 2008] |Bovy et al.||2012a|b|c| ), or if the structure 
varies with lo cation in the Galaxy. 

Meanwhile, |Nidever et al.|p014) found that the posi¬ 
tion of the locus of low-[a/Fe] (thin di sk) stars depends 
on location within the Galaxy (see also Edvardsson et al. 
1993), and it is possible that in the inner Galaxy the 
high- and low-[a/Fe] populations are connected, rather 
than distinct. Most previous observations were confined 
to the solar neighborhood, and use height above the 
plane or kinematics to separate between thick and thin 
disk populations. However , kinematical selections often 
misidentify stars (Bensby et al. 2014), and can remove 
intermediate or transitional populations; which may bias 
results flBovy et aL||2012c ). 

Observations of the metallicity distribution function 
(MDF) at different locations in the Galaxy can provide 
information about the evolutionary history across the 
disk. The MDF has generally only been well charac¬ 
terized in the solar neighbo rhood (e.g., |van den B ergh 


1962; Nordstr om et al.||2004| |Ak et al.||20071 ICa sagrande 

et al.||20fT| ~ Siebert et al. 20111 Schlesinge r et al.| 20 T 2 F 


_ ager et _ 

and in tKeTGalactic bulge"] Zoccali et al. 2008; Gonzalez 

et al. 2013; N ess et al.|2013 ) . The first observations of the 


MDF outside of the solar neighborhood were made us¬ 
ing APOGEE observations, and found differen ces i n Hie 


MDF as a function of Galactocentric radius (Anders et al 


2014). Metallicity distribution functions have long been 


used to constrain models of che mical evolution. Early 
chemical evolu tion models (e.g., Schmidt || 1963 Pagel & 


Patchett |1975| ) that attempted to explain the observed 


metal distribution in local G-dwarfs were simple, closed- 
box systems (no gas inflow or outflow) that over pre¬ 
dicted the number of metal-poor stars relative to obser¬ 
vations. This result commonly known as the “G-dwarf 


problem” 

(e.g., Pagel & Patchett ( 

1975); Rocha-Pinto 

& Maciel 

(1996|; |Schlesinger et al. 

|2012 

)). Solutions 


(e.g., Pagel 1997); observations of the MDF led to the 
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realization that gas dynamics play an important role in 
the chemical evolution of galaxies. However, it is not 
clear if the G-dwarf problem exists at all locations in the 
Galaxy, as there have been limited observations outside 
of the solar circle. 

Simulations and models of the chemical and kinemati- 
cal evolution of the Milky Way have become increasingly 


sophisticated (e.g., Hou et al. 2000; Chiap] 

pini et al.|2001; 

Schonrich & Binney||2009; Kubryk et al. 

2013 ; |Minchev 

et al. 
and o 

2013) and attempt to explain bot 
[ynamical history of the Galaxy. R 

i the chemical 
ecent chemical 


evolu tion models (e.g., Hon et ah||2000 Chiappini et al. 


2001) treat the chemistry of gas and stars in multiple el- 


ements across the entire disk, rather than just the solar 
neighborhood. Se veral s imu lations an d models find t hat 
“insid e-out” (e.g., Larson 1976 Kobayashi & Nakasato 
2011 [) and ^ “upsi de-down” (e.g., |Bournaud et al" 


2UU9 


2013) formation of the Galactic disk repro 


Bird et al 

duces observed trends in the Galaxy such as the radial 
gradient and the lower vertical scale heights of progres¬ 
sively younger populations. Centrally concentrated hot 
old disks, as seen in cosmological simulations, would re¬ 
sult in a decrease in scale-height with radius, which is 
not observed. In order to explain the presence of stars at 
high altitudes in the outer disk, in an inside-out forma¬ 
tion scenario, Minchev et al. (2015) suggested that disk 
flaring of mono-age populations is responsible. Such a 
view also explains the inversion of metallicity and [al- 
pha/Fe] gradients when the ve rtical distance from the 
disk midplane is increased (e.g ., Boeche et al.||2013| |An-| 
|ders et al.||2014| |Hayden et al.||2014|. Alternatively, the 
larger scale heights of older populations could be a con¬ 
sequence of satellite mergers. 

The radial mixing of gas and stars from their original 
birth radii has also been proposed as an important pro 


len et al. 1996 

Sellwood & Binney 2002 Roskar et al. 

2008; Schonric 

i & Binney| 2009; 

Loebman et al. 2011| 


curs through blurring, in which stars have increasingly 
more eccentric orbits and therefore variable orbital radii, 
and churning, where stars experience a change in angu¬ 
lar momentum and migrate to new locations while main¬ 
taining a circular orbit. However, there is much debate 
on the relative strength of mixing processes throughout 
the disk. Recent observations and modeling of the so¬ 
lar neighborhood have suggested that the local chemical 
structure of the dis k can be explained by blurring alone 
( Snaith et al. 2Q 14|, and tha t churning is not required, 


but see |Minchev et al. (2014). 


To disentangle these multiple processes and character¬ 
ize the history of the Milky Way disk, it is crucial to map 
the distribution of elements throughout the disk, beyond 
the solar neighborhood. This is one of the primary goals 
of the SDSS-III/APOGEE survey observed 146,000 stars 
across the Milky Way during three years of operation. 
APOGEE is a high-resolution (R^22,500) spectrograph 
operating in the iL-band, where extinction is 1/6 that 
of the V band. This allows observations of stars lying 
directly in the plane of the Galaxy, giving an unprece¬ 
dented coverage of the Milky Way disk. The main sur¬ 
vey goals were to obtain a uniform sample of giant stars 
across the disk with high resolution spectroscopy to study 
the chemical and kinematical structure of the Galaxy, in 


particular the inner Galaxy where optical surveys cannot 
observe efficiently due to high extinction. The APO GEE 
sur vey pr ovide s an RV precision of ^ 100 m s -1 (Nide- 
ver et al.|2015 ), and chemical abundances to within 0.1- 


0.2 dex for 15 different chemical elements (Garcia Perez 


et al. 


2015), in addition to excellent spatial coverage of 


the Milky Way from the bulge to the edge of the disk. 

In this paper we present results from the Twelfth Data 
Release (DR 12; Alam et al.|2015| of SDSS-III/APOGEE 
on the distribution!)! stars in the [a/Fe] vs. [Fe/H] plane 
and their metallicity distribution functions, across the 
Milky Way and at a range of heights above the plane. 
In Section 2 we discuss the APOGEE observations, data 
processing, and sample selection criteria. In Section 3 we 
present our observed results for the distribution of stars 
in [a/Fe] vs. [Fe/H] plane and metallicity distribution 
functions. In Section 4 we discuss our findings in the 
context of chemical evolution models. In the Appendix 
we discuss corrections for biases due to survey targeting, 
sample selection, population effects, and errors in the 
[a/Fe] determination. 

2. DATA AND SAMPLE SELECTION 

Data are taken from DR12, which contains stellar spec¬ 
tra and derived stellar parameters for stars observed dur¬ 
ing the three years of AP OGEE. APOGEE is on e of the 
main SDSS-III surveys (Eisenstein et al.||2011|), which 


uses the SDSS 2.5m telescope (Gunn et a. [.12DDB) i to ob¬ 
tain spectra for hundreds of stars per exposure? These 
stars cover a wide spatial extent of the Galaxy, and span 
a range of magnitudes between 8 < H < 13.8 for pri¬ 
mary science targets. Target selection is descr ibed in 
detai l in the APOGEE targeting paper (Zasowski et al . 
2013) and the APOGEE DR10 paper ( Ahn et al.||2014| ) 7 


Extinction and dereddening for each individual star is de¬ 
termine d using the Rayleigh-J eans Color Excess method 
(RJCE, Majewski et al. 2011), which uses 2MASS pho¬ 
tometry^ ([3km^kKn3On5006 ) in conjunct ion with near- 


IR photometry from the Spitzer/IRAC dFazio fc Team 


20041) GLIMP SE surveys (Benjamin et al.T2003; Ch urch- 
well et a~ [2009 ) where available, or from WISE ( Wright 
et al.||201(j|. In-depth discussio n of observing and reduc 


tion procedur es is described in (Hayden et al.||2014), the 
DR10 paper (A hn et al.||2014|), the APOGEE reduction 
pipeline pa per ( | N idever et al. 2015]), the DR12 calibra¬ 
tion p aper ([H oltzman et al. 2015|), the APOGEE linelist 
paper (| Zamora et al.|2015p, and the APOGEE Stellar Pa¬ 
rameters and Chemical A bundances Pipeline (ASPCAP, 
Garcia Perez et al.||2015) paper. 

For this paper, we select cool (T e ff < 5500 K) main 
survey (e.g., no ancillary program or Kepler field) giant 
stars (1.0 < log# < 3.8) with S/N> 80. Additionally, 
stars flagged as “Bad” due to being near the spectral 
library grid edge(s) or having poor spectral fits are re¬ 
moved. The cuts applied to the H-R diagram for DR12 
are shown in Figure [l] ASPCAP currently has a cutoff 
tem perature of 3500K on the cool side of the spect ral grid 
(see jGarcia Perez et al.|2015 Zamora et al.|2015 ), which 
could potentially bias our results against metal-rich stars. 
We correct for this metallicity bias by imposing the lower 
limit on surface gravity of logg > 1.0, as this mitigates 
much of the potential bias due to the temperature grid 
edge in the observed metallicities across the disk; for a 
detailed discussion see the Appendix. Applying these re- 
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strictions to the DR12 catalog, we have a main sample of 
69,919 giants within 2 kpc of the midplane of the Milky 
Way. We restrict our study to the disk of the Milky Way 
(R > 3 kpc). For a detailed discussion of the M DF of 
the Galactic bu l ge with APOGEE observations see | Gar- 1 


cia Perez et al. 


(in prep). 

[Fe/HJ has been calibrated to literature values for a 
large set of reference stars and clusters, and all other 
spectroscopic parameters are calibrated using clus ters to 
remove trends w ith temperature, as described in |Holtz- 
man et al. (2015). These corrections are similar to those 


applied to the DH10 sample by Meszaros et al. (2013), as 
there are slight systematic offsets observed in the ASP- 
CAP parameters compared to reference values, and in 
some cases trends with other parameters (e.g., abun¬ 
dance trends with effective temperature). The accuracy 
of the abundances has improved from DR10, in particular 
for [a/Fe] , as self-consistent model atmospheres rather 
than scaled solar atmospheres were used in the latest 
data release (see [Zamora et al.||2015|), improving the ac¬ 
curacy of many parameters. The typical uncertainties in 
the spectroscopic parameters from||Holtzman et al. ([2015 ) 
are: 0.11 dex in logg, 92 K in T e ^, and 0.05 dex in [Fe/HJ 
and [a/Fe] . 

Giants are not perfectly representative of underlying 
stellar populations, as they are evolved stars. There is a 
bias against the oldest populations when using giants as a 
tracer population, with the relative po pulation sampling 
being oc r -0,6 , where r is age (e.g., (Girardi & Salaris 
2001|)). It is difficult to correct for this population sam¬ 
pling effect, as it depends on the detailed star formation 
histories at different locations throughout the Galaxy, 
and likely requires detailed population synthesis models. 
Because of this, we do not correct for the non-uniform 
age sampling of giants and our MDFs are slightly biased 
against the oldest (and potential more metal-poor) stars 
of the disk. For additional discussion on population sam¬ 
pling, see the Appendix. 


2.1. Distances 

Distances for each star are determined from the de¬ 
rived stellar parameters an d PARSEC isochrone s from 
the Padova-Trieste group (Bressan et al. 2012) based 
on Bayesian statis t ics, following methods des cribe d by 
Burnett & Binneyl (|2010|)|Burnett et al.| (|2011|),and|Bin- 
ney et al. (2011); see also Santiago et al. r Pl5l )- TEe 
isochrones range in metallicity from —2.5 < [Fe/H] < 
+0.6, with a spacing of 0.1 dex, and ages (r) ranging 
from 100 Myr to 20 Gyr with spacing of 0.05 dex in log r. 
We calculate the probability of all possible distances us¬ 
ing the extinction-corrected magnitude (from RJCE, as 
referenced above), the stellar parameters from ASP CAP, 
and the PARSEC isochrones using Bayes’ theorem: 


P (model\data) = 


P (data | model ) P (model ) 
P(data) 


where model refers to the isochrone parameters (T e ff, 
logg, r, [Fe/H] , etc.) and physical location (in our case 
the distance modulus). Data refers to the observed spec¬ 
troscopic and photometric parameters for the star. For 
our purposes, we are interested in the distance modulus 
(/i) only, so P(model|data) is: 



- 1.0 - 0.5 0.0 0.5 

< [Fe/H] > 


Figure 1 . The spectroscopic H-R diagram for the full calibrated 
APOGEE sample, where the mean metallicity in each log g,T e ff 
bin is shown. The gray box denotes the selected sample of 69,919 
stars presented in this paper. 


P(model\data) = P{/i) oc 


n 


exp 



dl 


where o j is the observed spectroscopic parameter, Ij is 
the corresponding isochrone parameter, and a Q . is the er¬ 
ror in the observed spectroscopic parameter. Additional 
terms can be added if density priors are included, but we 
did not include density priors for the distances used in 
this paper; our effective prior is flat in distance modulus. 
The distance modulus most likely to be correct given the 
observed parameters and the stellar models is determined 
by creating a probability distribution function (PDF) of 
all distance moduli. We use isochrone points within 3a 
of our observed spectroscopic temperature, gravity, and 
metallicity to compute the distance moduli, where the 
errors in the observed parameters are given in the data 
section above. To generate the PDF, the equation above 
is integrated over all possible distance moduli, although 
in practice we use a range of distance moduli between the 
minimum and maximum magnitudes from the isochrone 
grid matches to reduce the required computing time. 

The peak, median, or average of the PDF can be used 
to estimate the most likely distance modulus for a given 
star. For this paper, we use the median of the PDF 
to characterize the distance modulus. The error in the 
distance modulus is given by the variance of the PDF: 




< n > 2 


The radial distance from the Galactic center is computed 
assuming a solar distance of 8 kpc from the Galactic 
center. Distance accuracy was tested by comparing to 
clusters observed by APOGE E and using simulate d ob¬ 
servations from TRILEGAL (Girardi et al. 2005). On 
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4 6 8 10 12 14 

R(kpc) 


Figure 2. The Galactic R-z distribution for the sample of 69,919 
stars used in this analysis. R is the projected planar distance 
from the Galactic center, while z is the distance from plane of the 
Galaxy. Each star is plotted at the location implied by the median 
of its distance modulus PDF. 


average, the distances are accurate at the 15 — 20% level. 
A more detailed discussio n of the distances can be found 
in|Holtzman et al. (2015). The stellar distribution in the 
R-z plane for the sample used in this paper is shown in 
Figure |2j 


3. RESULTS 
3.1. [a/Fe] vs. [Fe/H] 

We present results of the distribution of stars in the 
[a/Fe] vs. [Fe/H] plane in the solar neighborhood (7 < 
R < 9 kpc, 0 < \z\ < 0.5 kpc) in Figure |3l The stel¬ 
lar distribution in the [a/Fe] vs. [Fe/H] plane in the 
solar neighborhood is characterized by two distinct se¬ 
quences, one starting at high-[a/Fe] and the other at ap¬ 
proximately solar-[a/Fe] abundances. We use the term 
sequences or tracks to describe the behavior of the low 
and high-[a/Fe] populations, but this description does 
not necessarily imply the sequences to be evolutionary in 
nature. The high-[a/Fe] sequence has a negative slope, 
with the [a/Fe] ratio decreasing as [Fe/H] increases, and 
eventually merges with the low-[a/Fe] sequence at [Fe/H] 
~ +0.2. The low-[a/Fe] sequence has a slight decrease 
in [a/Fe] abundance as metallicity increases, except for 
the most metal-rich stars ([Fe/H] >0.2) where there the 
trend flattens. The lower envelope of the distribution has 
a concave-upward, bowl shape. However, these trends 
are small, and the sequence is within ~ 0.1 dex in a 
abundance across nearly a decade in metallicity. It is 
unclear which sequence the most metal-rich stars belong 
to in the solar neighborhood, because the low- and high- 
[a/Fe] sequences appear to merge at these super-solar 


studies of the solar neighborhood (e.g., Adibekyan et al. 

2012; IRamfrez et al. 20131 IBensby et al.| 

20l4| INidever 

et al. 2014 Recio-Blanco et al. 2014]), whic 

i also find two 
plane. While 

distinct sequences in the [a/Fe J vs. [Fe/H 


some surveys targeted kinematic “thin” and “thick” disk 
samples in a wa y that could amplify bim odality, our sam¬ 
ple (and that of Adibekyan et al.|2011 ) has no kinematic 


selection. 

APOGEE allows us to extend observations of the dis¬ 


tribution of stars in the [a/Fe] vs. 
much of the disk (3 < R < 15 kpc, \z 


Fe/H] plane across 
< 2 kpc) as shown 


in Figure [4|. The most striking feature of the stellar 
distribution in the [a/Fe] vs. [Fe/H] plane in the inner 
disk (3 < R < 5 kpc) is that the separate low-[a/Fe] 
sequence evident in the solar neighborhood is absent— 
there appears to be a single sequence starting at low 
metallicities and high-[a/Fe] abundances, which ends at 
approximately solar-[a/Fe] and high metallicity ([Fe/H] 
~ +0.5). The metal-rich solar-[a/Fe] stars dominate the 
overall number density of stars close to the plane. These 
stars are confined to the midplane, while for \z\ > 1 kpc 
the majority of the stars in the inner Galaxy are metal- 
poor with high-[a/Fe] abundances. 

In the 5 < R < 7 kpc annulus, the population of low- 
[a/Fe] , sub-solar metallicity stars becomes more promi¬ 
nent, revealing the two-sequence structure found in the 
solar neighborhood. However, the locus of the low alpha 
sequence is significantly more metal rich ([Fe/H]~ 0.35) 
than our local sample. The mean metallicity of the 
low alpha sequence and its dependence on radius largely 
drives the observed Galactic metallicity gradients. As 
|z| increases, the relative fraction of high- and low-[a/Fe] 
stars changes; for |z| > 1 kpc, the bulk of the stars belong 
to the high-[a/Fe] sequence and have sub-solar metallic- 
ities. 

Towards the outer disk (R > 9 kpc), the locus of the 
low- [a/Fe] sequence shifts towards lower metallicities. 
Much like the rest of the Galaxy, the low-[a/Fe] sequence 
dominates the number density close to the plane of the 
disk. As \z\ increases, the high-[a/Fe] fraction increases 
(for 9 < R < 13 kpc), but it never becomes the dominant 
population at high \z\ as it does in the solar neighbor¬ 
hood or inner disk. For R > 13 kpc, there are almost no 
high-[a/Fe] stars present; at all heights above the plane 
most stars belong to the low-[a/Fe] sequence. For R > 11 
kpc, the relative number of super-solar metallicity stars 
is low compared to the rest of the Galaxy; these stars are 
confined to the inner regions (R < 11 kpc) of the disk. 
The spread in metallicity for the very outer disk (R > 13 
kpc) is small: most stars are within [Fe/H] ^ —0.4 ±0.2 
d ex at all heig hts about the plane. 

INideve r et al.| (|2014|) used the APOGEE Red Clump 
Catalog (Hovy et al. "[2014 ) to categorize the distribu¬ 
tion of stars in the [a/Fe] v s. [Fe/H] plane outs ide of 
the solar neighborhood. The Nidever et al. (2014|) sam¬ 
ple is a subset of the same data presented in this paper. 
The red clump offers more precise distance and abun¬ 
dance determinations compared to the entire DR12 sam¬ 
ple, but it covers a m ore restric ted distance and metallic- 
ity range. Nidever et al. (2014) find that the high-[a/Fe] 
sequence found in the solar neighborhood is similar in 
shape in all areas of the Galaxy where it could be ob¬ 
served (5 < R < 11 kpc, 0 < \z\ < 2 kpc). Here we 
expand these observations to larger distances and find a 
similar result; the shape of the high-[a/Fe] sequence does 
not vary significantly with radius, although the very in¬ 
ner Galaxy does show a hint of small differences. There 
appears to be a slight shift towards lower-[a/Fe] for the 
same metallicities by 0.05 dex compared to the high- 
[a/Fe] sequence observed in the rest of the disk. This 
variation may be caused by temperature effects; the stars 
in the inner Galaxy are all cool (T e ff < 4300 K), and 
there is a slight temperature dependence of [a/Fe] abun¬ 
dance for cooler stars, as discussed further in the Ap¬ 
pendix. Although the high-[a/Fe] sequence appears sim- 
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Solar Neighborhood 



[Fe/H] 


Figure 3. The observed [a/Fe] vs. [Fe/H] distribution for the solar neighborhood (7 < R < 9 kpc, 0 < |^| < 0.5 kpc). The left panel is 
down-sampled and shows only 20% of the observed data points in the solar circle. The right panel shows the entire sample in the solar 
neighborhood, with contours denoting 1, 2, 3cr of the overall densities. There are two sequences in the distribution of stars in the [a/Fe] vs. 
[Fe/H] plane, one at solar-[a/Fe] abundances, and one at high-[a/Fe] abundances that eventually merges with the solar-[a/Fe] sequence at 
[Fe/H] - 0.2. 



[Fe/H] 


Figure 4. The stellar distribution of stars in the [a/Fe] vs. [Fe/H] plane as a function of R and \z\. Top: The observed [a/Fe] vs. [Fe/H] 
distribution for stars with 1.0 < \z\ < 2.0 kpc. Middle: The observed [a/Fe] vs. [Fe/H] distribution for stars with 0.5 < \z\ < 1.0 kpc. 
Bottom: The observed [a/Fe] vs. [Fe/H] distribution for stars with 0.0 < \z\ < 0.5 kpc. The grey line on each panel is the same, showing 
the similarity of the shape of the high-[a/Fe] sequence with R. The extended solar-[a/Fe] sequence observed in the solar neighborhood is 
not present in the inner disk (R < 5 kpc), where a single sequence starting at high-[a/Fe] and low metallicity and ending at solar-[a/Fe] 
and high metallicity fits our observations. In the outer disk (R > 11 kpc), there are very few high-[a/Fe] stars. 


ilar at all observed locations, as noted above the number 
of stars along the high-[a/Fe] sequence begins to decrease 
dramatically for R > 11 kpc: there are almost no stars 
along the high-[a/Fe] sequence in the very outer disk 
(13 < R < 15 kpc). 

To summarize our results for the distribution of stars 
in the [a/Fe] vs. [Fe/H] plane: 

• There are two distinct sequences in the solar neigh¬ 
borhood, one at high-[a/Fe] , and one at solar- 
[a/Fe] , which appear to merge at [Fe/H] ~ +0.2. 
At sub-solar metallicities there is a distinct gap be¬ 
tween these two sequences. 

• The abundance pattern for the inner Galaxy can 
be described as a single sequence, starting at low- 
metallicity and high-[a/Fe] and ending at approxi¬ 
mately solar-[a/Fe] and [Fe/H] = +0.5. The most 
metal-rich stars are confined to the midplane. 


• The high-[a/Fe] sequence appears similar at all lo¬ 
cations in the Galaxy where it is observed (3 < 
R < 13 kpc). 

• Stars with high-[a/Fe] ratios and the most metal- 
rich stars ([Fe/H] >0.2) have spatial densities that 
are qualitatively consistent with short radial scale- 
lengths or a truncation at larger radii and have low 
number density in the outer disk. 

• The relative fraction of stars between the low- and 
high-[a/Fe] sequences varies with disk height and 
radius. 

3.2. Metallicity Distribution Functions 

With three years of observations, there are sufficient 
numbers of stars in each Galactic zone to measure MDFs 
in a number of radial bins and at different heights above 
the plane. We present the MDFs in radial bins of 2 kpc 


Density Percentile 
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Table 1 

Metallicity Distribution Functions in the Milky Way 


R Range (kpc) 

N* 

<[Fe/H] > Peak [Fe/H] 

^[Fe/H] 

Skewness 

Kurtosis 

1.00 < \z\ < 2.00 

3 < R < 5 

465 

-0.42 

-0.27 

0.29 

-0.48+0.14 

1.05+0.28 

5 < R < 7 

846 

-0.36 

-0.33 

0.29 

-0.32+0.13 

1.27+0.26 

7 < R < 9 

4136 

-0.31 

-0.27 

0.28 

-0.53+0.06 

1.52+0.17 

9 < R < 11 

1387 

-0.29 

-0.27 

0.25 

-0.37+0.13 

1.69+0.36 

11 < R < 13 

827 

-0.29 

-0.38 

0.23 

-0.40+0.21 

2.58+0.62 

13 < R < 15 

207 

-0.39 

-0.43 

0.17 

-0.60+0.73 

3.84+3.16 

0.50 < \z\ < 1.00 

3 < R < 5 

841 

-0.19 

-0.33 

0.32 

-0.50+0.11 

0.50+0.31 

5 < R < 7 

1408 

-0.12 

-0.18 

0.29 

-0.50+0.09 

0.39+0.35 

7 < R < 9 

4997 

-0.10 

-0.02 

0.25 

-0.49+0.06 

0.67+0.22 

9 < R < 11 

3702 

-0.15 

-0.23 

0.21 

-0.22+0.10 

0.99+0.47 

11 < R < 13 

2169 

-0.23 

-0.27 

0.19 

+0.28+0.11 

0.92+0.39 

13 < R < 15 

568 

-0.33 

-0.33 

0.19 

-0.60+0.39 

4.63+1.50 

0.00 < \z\ < 0.50 

3 < R < 5 

2410 

+0.08 

+0.23 

0.24 

-1.68+0.12 

4.01+0.83 

5 < R < 7 

5195 

+0.11 

+0.23 

0.22 

-1.26+0.08 

2.53+0.52 

7 < R < 9 

13106 

+0.01 

+0.02 

0.20 

-0.53+0.04 

0.86+0.26 

9 < R < 11 

19930 

-0.11 

-0.12 

0.19 

-0.02+0.03 

0.49+0.14 

11 < R < 13 

6730 

-0.21 

-0.23 

0.18 

+0.17+0.06 

0.79+0.21 

13 < R < 15 

912 

-0.31 

-0.43 

0.18 

+0.47+0.13 

1.00+0.29 


Note. — Statistics for the different MDFs in the across the Milky Way disk. The kurtosis 
is defined as the fourth standardized moment-3, such that a normal distribution has a kurtosis 
of 0. 


Table 2 

a Distribution Functions in the Milky Way 


R Range (kpc) 

N* 

< [a/ H]> Peak [a/ H] 

°Wh] 

Skewness 

Kurtosis 

1.00 < |z| < 2.00 

3 < R < 5 

468 

-0.20 

-0.08 

0.26 

-1.12+0.23 

3.02+0.91 

5 < R < 7 

853 

-0.16 

-0.08 

0.27 

-1.25+0.18 

4.17+0.67 

7 < R < 9 

4145 

-0.14 

-0.08 

0.24 

-1.11+0.09 

3.70+0.34 

9 < R < 11 

1393 

-0.16 

-0.18 

0.23 

-0.93+0.22 

4.68+0.90 

11 < R < 13 

831 

-0.20 

-0.27 

0.22 

-0.82+0.29 

4.61+1.06 

13 < R < 15 

210 

-0.30 

-0.38 

0.21 

-1.91+0.57 

9.19+1.80 

0.50 < \z\ < 1.00 

3 < R < 5 

844 

-0.03 

+0.02 

0.27 

-1.12+0.20 

3.34+0.87 

5 < R < 7 

1409 

+0.02 

+0.07 

0.23 

-0.73+0.15 

1.73+0.74 

7 < R < 9 

4997 

-0.01 

+0.02 

0.21 

-0.39+0.06 

0.78+0.25 

9 < R < 11 

3703 

-0.08 

-0.12 

0.19 

-0.05+0.10 

0.98+0.48 

11 < R < 13 

2170 

-0.16 

-0.23 

0.18 

0.34+0.19 

1.74+1.00 

13 < R < 15 

568 

-0.25 

-0.27 

0.17 

-0.42+0.45 

4.65+1.88 

0.00 < \z\ < 0.50 

3 < R < 5 

2414 

+0.16 

+0.27 

0.20 

-1.91+0.20 

7.17+1.53 

5 < R < 7 

5195 

+0.16 

+0.23 

0.19 

-1.07+0.10 

2.54+0.78 

7 < R < 9 

13109 

+0.06 

+0.07 

0.18 

-0.34+0.07 

1.20+0.50 

9 < R < 11 

19930 

-0.07 

-0.08 

0.17 

+0.27+0.03 

0.35+0.11 

11 < R < 13 

6731 

-0.16 

-0.23 

0.16 

+0.38+0.07 

1.17+0.37 

13 < R < 15 

914 

-0.25 

-0.33 

0.17 

+0.04+0.39 

3.75+1.82 


Note. — Statistics for the different ADFs in the across the Milky Way disk. The 
kurtosis is defined as the fourth standardized moment-3, such that a normal distribution 
has a kurtosis of 0. 


between 3 < R < 15 kpc, and at a range of heights above 
the plane between 0 < |z| < 2 kpc, in Figure [5] The 
MDFs are computed with bins of 0.05 dex in [Fe/H] for 
each zone. Splitting the sample into vertical and radial 
bins allows us to analyze the changes in the MDF across 
the Galaxy, but also minimizes selection effects due to 
the volume sampling of the APOGEE lines of sight and 
our target selection. 

Close to the plane (top panel of Figure H hi < °- 5 
kpc) radial gradients are evident throughout the disk. 
The peak of the MDF is centered at high metallicities 
in the inner Galaxy ([Fe/H] = 0.32 for 3 < R < 5 kpc), 
roughly solar in the solar neighborhood (M/H]=+0.02 


for 7 < R < 9 kpc), and low metallicities in the outer 
disk ([Fe/H] = —0.48 for 13 < R < 15 kpc). The radial 
gradients observed in the MDF are similar to those mea¬ 
sured across the disk in previous studies with APOGEE, 


as is the shift in the peak of the MDF (e.g., Anders et al. 
2014 Hayden et al.||2014 ) 


The most striking feature of the MDF close to the plane 
is the change in shape with radius. The inner disk has a 
large negative skewness (—1.68 ±0.12 for 3 < R < 5 kpc, 
see Table [l]), with a tail towards low metallicities, while 
the solar neighborhood is more Gaussian in shape with 
a slight negative skewness (—0.53 ± 0.04), and the outer 
disk is positively skewed with a tail towards high metal- 

























Figure 5. The observed MDF for the entire sample as a function 
of Galactocentric radius over a range of distances from the plane. 
The shape and skewness is a function of radius and height. Close 
to the plane, the inner Galaxy (3 < R < 5 kpc) is a negatively 
skewed distribution and a peak metallicity at ~ 0.25 dex, while 
the outer disk (R > 11 kpc) has a positively skewed distribution 
with peak metallicity of ~ —0.4 dex . For \z\ > 1 kpc, the MDF is 
fairly uniform at across all radii. 

licities (+0.47T0.13 for 13 < R < 15 kpc). The shape of 
the observed MDF of the solar neighborhood is in goo d 
agreement with the MDF measured by the GCS ([Nord¬ 

stro m et al.|2004[ |Holmberg et al.|2 007 ; Casagrande et al. 
201 ip, who measure a peak of just below solar metallicity 
anclalso find a similar negative skewness to our obser¬ 
vations. There is a slight offset in the peak metallicity, 
with the APOGEE observations being more metal-rich 
by ~ 0.1 dex, but the shapes are extremely similar. Close 
to the plane, the distributions are all leptokurtic, with 
the inner Galaxy (3 < R < 7 kpc) being more strongly 
peaked than the rest of the disk. 

As \z\ increases, the MDF exhibits less variation with 
radius. For \z\ > 1 kpc (bottom panel of Figure [5]), the 
MDF is uniform with a roughly Gaussian shape across 
all radii, although it is more strongly peaked for the very 
outer disk (R > 13 kpc). However, the populations com¬ 
prising the MDF(s) at these heights are not the same. In 
the inner disk, at large heights above the plane the high- 
[a/Fe] sequence dominates the number density of stars. 
In the outer disk (R> 11 kpc), the stars above the plane 
are predominantly solar-[a/Fe] abundance. The unifor¬ 
mity of the MDF at these large heights is surprising given 
the systematic change in the [a/Fe] of the populations 
contributing to the MDF. The MDF is similar for all 
heights above the plane in the outer disk (R > 11 kpc). 
At these larger heights above the plane, the MDFs are 
leptokurtic as well, but the trend with radius is reversed 
compared to the distributions close to the plane. For 
\z\ > 1 kpc, the distributions in the outer disk (R > 11 



- 1.5 - 1.0 - 0.5 0.0 0.5 

[Fe/H] 


Figure 6. The observed MDF for stars with [a/Fe] > 0.18 as a 
function of Galactocentric radius over a range of distances from 
the plane. There is little variation in the MDF with Galactocentric 
radius. There is a shallow negative vertical gradient, as stars with 
\z\ >1 kpc are slightly more metal-poor than stars close to the 
plane. 

kpc) are more strongly peaked than the MDFs for the 
rest of the disk. 

As noted above, the high-[a/Fe] sequence is fairly con¬ 
stant in shape with radius. The MDF for stars with high- 
[a/Fe] abundance ([a/Fe] > 0.18) is presented in Figure 
p] The high-[a/Fe] sequence appears uniform across the 
radial range in which it is observed (3 < R < 11 kpc); but 
does display variation with height above the plane. Close 
to the plane, the MDF is peaked at [Fe/H] = —0.3 over 
the entire radial extent where there are large numbers of 
high- [a/Fe] stars; and the shape is also the same at all 
locations. However, as \z\ increases, the MDF shifts to 
slightly lower metallicities, with a peak [Fe/H] = —0.45 
for \z\ > 1 kpc. At any given height, there is little vari¬ 
ation in the shape or peak of the MDF with radius for 
these high-[a/Fe] stars. 

Simple chemical evolution models often use instanta¬ 
neous recycling approximations where metals are imme¬ 
diately returned to the gas reservoir after star formation 
occurs. This approach may not be a good approxima¬ 
tion for all chemical elements in the stellar population, 
but it is more accurate for a elements, which are pro¬ 
duced primarily in SNII. We present the distribution of 
[a/H] (ADF) across the disk in Figure [ 7 ] The ADF is 
similar in appearance to the MDF, with radial gradients 
across the disk and the most a-rich stars belonging to 
the inner Galaxy. The observed change in skewness with 
radius of the ADF (Table[2| is similar to that of the MDF 
(Table [l). Close to the plane, the ADF of the inner disk 
(3 < R < 5 kpc) is more strongly peaked than the MDF 
in the same zone, with a significantly larger kurtosis. The 
main difference between the ADF and the MDF is above 
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[ct/H] 


Figure 7. The observed distribution of [ck/H] for the entire 
sample as a function of Galactocentric radius over a range of 
distances from the plane. The results are quite similar to that 
of the MDF, except for stars with \z\ > 1 kpc and R < 9 kpc, 
where the ADF has larger abundances than the MDF at the same 
locations. 


the plane of the disk, we observe radial gradients in the 
ADF at all heights above the plane, which is not the case 
for the MDF. At large \z\, the ADF is significantly more 
positive than the MDF from the inner Galaxy out to the 
solar neighborhood by ~ 0.25 dex. This result does not 
hold true in the outer disk above the plane (R > 11 kpc, 
\z\ > 1 kpc), where the ADF again has similar abundance 
trends to that of the MDF. 

To summarize our results for the MDFs across the disk: 

• Metallicity gradients are clearly evident in the 
MDFs, with the most metal-rich populations in the 
inner Galaxy. 

• The shape and skewness of the MDF in the mid¬ 
plane is strongly dependent on location in the 
Galaxy: the inner disk has a large negative skew¬ 
ness, the solar neighborhood MDF is roughly Gaus¬ 
sian, and the outer disk has a positive skewness. 

• The MDF becomes more uniform with height. For 
stars with \z\ > 1 kpc it is roughly Gaussian with a 
peak metallicity of [Fe/H] ~ —0.4 across the entire 
radial range covered by this study (3 < R < 15 
kpc). 

• The MDF for the outer disk (R > 11 kpc)is uni¬ 
form at all heights above the plane, showing less 
variation in metallicity with \z\ than the rest of the 
disk. 

• The MDF for stars with [a/Fe] > 0.18 is uniform 
with i?, but has a slight negative vertical gradient. 


• The ADF has many of the same features observed 
in the MDF, but shows differences for stars out of 
the plane (\z\ > 1 kpc) and with R < 9 kpc, where 
stars tend to have higher [cu/H] than [Fe/H] . 


4. DISCUSSION 

Consistent with previous studies, we find that: (a) the 
solar-neighborhood MDF is approximately Gaussian in 
[Fe/H], with a peak near solar metallicity, (b) the dis¬ 
tribution of stars in the [a/Fe] -[Fe/H] plane is bimodal, 
with a high-[a/Fe] sequence and a low-[a/Fe] sequence, 
(c) the fraction of stars with high [a/Fe] increases with 
| z |, and (d) there is a radial gradient in the mean or me- 
d ian value of [Fe/H] for stars near the midplane. Like 
Ni dever et al.| (2014), we find that the location of the 
high-[a/Fe] sequence is nearly independent of radius and 
height above the plane, a result we are able to extend 
to larger and smaller R and to lower metallicity. Two 
striking new results of this study are (e) that the [a/Fe] 
-[Fe/H] distribution of the inner disk (3 < R < 5 kpc) is 
consistent with a single evolutionary track, terminating 
at [Fe/H] ~ +0.4 and roughly solar [a/Fe] , and (f) that 
the midplane MDF changes shape, from strong negative 
skewness at 3 < R < 7 kpc to strong positive skewness 
atll</£<15 kpc, with the solar annulus lying at the 
transition between these regimes. 

The midplane inner disk MDF has the characteristic 
shape predicted by one-zone chemical evolution models, 
with most stars formed after the ISM has been enriched 
to an “equilibrium” metallicity controlle d mai nly by Hie 
outflow mass loading parameter r] (see Andrews et al. 
2015 , for detailed discussion). Traditi onal closed box 
or leaky box models (see, e.g., §5.3 of Binney & Mer- 
rifield |1998]) are a limiting case of such models, with 
no accretion. These models predict a metallicity distri¬ 
bution dN/dlnZ oc Zexp(—Z/p e ff) where the effective 
yield is related to the IMF-averaged population yield by 
Peff = p/(l + Tf). The positively skewed MDFs of the 
outer disk could be a distinctive signature of radial mi¬ 
gration, with a high-metallicity tail populated by stars 
that were born in the inner Galaxy. The change of [a/Fe] 
distributions with height could be a consequence of heat¬ 
ing of the older stellar populations or of forming stars in 
progressively thinner, “cooler” populations as turbulence 
of the early sta r-forming disk decreases. As discussed by 
Nidever et al. 2014[), the constancy of the high- [a/Fe] 
sequence implies uniformity of the star formation effi¬ 
ciency and outflow mass loading during the formation of 
this population. 

In the next subsection we discuss the qualitative com¬ 
parison between our results and several recent models of 
Milky Way chemical evolution. We then turn to a more 
detailed discussion of radial migration with the aid of 
simple quantitative models. We conclude with a brief 
discussion of vertical evolution. 


4.1. Comparison to Chemical Evolution Models 

Metallicity distribution functions are useful observa¬ 
tional tools in constraining the chemical history of the 
Milky Way. The first chemical evolution models were 
simple closed-box systems, with no gas inflow or out¬ 
flow, and often employed approximations such as instan¬ 
taneous recycling. These models over-predicted the num¬ 
ber of metal-poor stars in the solar-neighborhood com- 
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pared to observations of G-dwarfs (e.g., Schmidt 1963 r 
Pagel fe Pat chet t 11975 ) ? a discrepancy known as the U G- 


Dwarf Problem”. These first observations made it clear 
that the chemical evolution of the solar neighborhood 
could not be described by a simple closed-box model; in¬ 
flow and outflow of gas, along with more realistic yields 
(i.e., no instantaneous recycling) were required to repro¬ 
duce observations. The MDF can therefore be used to 
inform and tune chemical evolution models and provide 
information such as the star formation history and rela¬ 
tive gas accretion or outflow rates at every location where 
the MDF can be measured. APOGEE observations pro¬ 
vide the first thorough characterizations of the MDF of 
the disk outside of the solar neighborhood, allowing a 
more complete characterization of the chemical evolu¬ 
tionary history of the Galaxy. 

Additions such as gas inflow and outflow to chemical 
evolution models have been able to better reproduce ob¬ 
servations of the solar neighborhood, in particular the 


MDF and stellar distribution [a/Fe] vs. 

'Fe/Hl i 

plane. 

The two-infall model from Chiappini et al. 

( 

1997 

2001) 

treats the disk as a series of annuli, into wJ 

mch gas ac- 


cretes. In this model, an initial gas reservoir forms the 
thick disk, following the high-[a/Fe] evolutionary track. 
As the gas reservoir becomes depleted, a lull in star for¬ 
mation occurs. SNela gradually lower the [a/Fe] ratio 
of the remaining gas reservoir as a second, more grad¬ 
ual infall of pristine gas dilutes the reservoir, lowering 
the overall metallicity but retaining the low-[a/Fe] abun¬ 
dance of the ISM. Once the surface density of the gas is 
high enough star formation resumes, forming the metal- 
poor end of the solar-[a/Fe] sequence. The MDF from 
the two-infall model is in general agreement with our 
observations of the so lar neighborhood (see Figure 7 of 
Chiappi ni et al.| 1997), with a peak metallicity near solar 
and a slight negative skewness towards lower metallici- 
ties. Additionally, this model reproduces general trends 
found in the distribution of stars in the [a/Fe] vs. [Fe/H] 
plane, in particular with the dilution of the metallicity 
of the existing gas reservoir with pristine gas to form the 
low-[a/Fe] sequence. This is one possible explanation for 
the observed low-[a/Fe] sequence: pristine gas accretes 
onto the disk and mixes with enriched gas from the inner 
Galaxy, keeping solar-[a/Fe] (or intermediate [a/Fe] that 
are later lowered to solar-[a/Fe] ratios by SNela) ratios 
but lowering the metallicity from +0.5 dex in the inner 
disk to the lower metallicities found in the outer Galaxy. 
This model, and its ability to explain phenomenologically 
both the MDF and stellar distribution in the [a/Fe] vs. 
[Fe/H] plane of the solar neighborhood, highlight the po¬ 
tential importance of gas flow in the evolutionary history 
of the Galaxy. 

T he che mical evolution model from ISchonrich ~fe Bin-1 
(2009) includes radial migration in the processes gov- 


ney 


erning"the evolution of the disk. Their models find that 
the peak of the MDF is a strong function of radius, with 
the inner Galaxy being more metal-rich than the outer 
Galaxy. The peaks of their MDFs are similar to those 
observed in APO GEE at different radii (see Figure 11 
of Schonrich & Binney 2009). However, their distribu¬ 
tions are significantly" more Gaussian and less skewed 
than we o bserve in our sample . The model distribu- 


ner Galaxy to the outer Galaxy, as we observe, but the 
magnitude of th e skewness is not large. The Scho nrich] 
& Binney ( |2009 ) results for [O/Fe] vs. [Fe/HJ is not in 


good agreement to the distribution of stars in the [a/Fe] 
vs. [Fe/H] plane ob served wit h the APOGEE sample 
(see Figure 9 of Schonrich & Binney] 2009). Their mod¬ 
els do not have two distinct sequences m this plane, but a 
more continuous distribution of populations that start at 
low-metallicity with high-[a/Fe] ratios and end at high- 
metallicity and solar-[a/Fe] , similar to many other mod¬ 
els with inside-out formation. Their model has have a 
much larger dynamic range than the APOGEE sample 
in [a/Fe] abundance, with their [O/Fe] ratio extending to 
~ +0.6, while the APOGEE [a/Fe] abundances extend 
only to ^ +0.3, with the thin disk sequence in APOGEE 
shifted to slightly higher metallicities and lower-[a/Fe] 
abundan ces than the thin-disk sequence presented by 


Schonrich & Binney (2009). 

Recent IN- body smooth particle hydrodynamics simu¬ 


lations from Kubryk et al. (2013) track the impact of 
migration on the stellar distrib ution of their disc ga laxy. 
Their MDFs (see Figure 10 oflKubryk et al.||2013) are 
fairly uniform throughout the disk, appearing similar to 
the MDF of the solar neighborhood in the APOGEE 
observations, with a peak n ear so lar abundances and a 
negative skewness, 
nificant shifts in t 


Kubryk et al. (2013) do not find sig- 
re peak or skewness with radius in 
their simulations, contrary to what is observed in the 
APOGEE observations. They postulate that the unifor¬ 
mity of the peak of the MDF is due to the lack of gas 
infall in the simulation, highlighting the importance that 
gas dynamics can play in the MDF across the disk. 

The most sophisticated chemical evolution models to 
date use cosmological simulations of a Milky Way ana¬ 
log Galaxy, and paint a chemical evolution model on top 
of the simul ated galaxy. Recent simulations by Mi nchev| 
et al. (2013) match many observed properties of the disk, 


such as the stellar distri bution in the [a/Fe] v s. [Fe/H] 
plane (see Figure 12 of|Minchev et al.| 2013|), and the 
flattening of the radial gradient with height (e.g., |An-| 
ders et aL][2 0 14l [Hayden et al.[|2014| see Figure 10 of 


Min chev et al.||2014 ). Their simulations also provide de 


tailed metallicity distributions for a large range of radii. 
We find that the MDFs from the simulations closely 
match the APOGEE observations of the MDF in the so¬ 
lar neighborhood and the u niformity of the MDF above 
the plane (see Figure 9 of Minchev et al.|2014). While we 
have g ood agreement with the MDFs lrom |lVlinchev et ah] 
(2014) above the plane at all radii, the MDFs from their 
simulation do not reproduce the change in the peak and 
skewness of the MDF observed close to the pla ne in the 
inner and o uter disk. The MDFs presented in Minchev] 


tions from Schonrich & Binney (2009) appear to have 
a shift from negative to positive skewness from the in- 


et al. (2014) have the same peak metallicity at all radii, 
which is not observed in APOGEE close to the plane. 
The metal-rich components of the MDFs from the sim¬ 
ulation are in the wings of the distributions, leading to 
positively skewed MDFs in the inner Galaxy, and roughly 
Gaussian shapes in the outer disk. The APOGEE ob¬ 
servations have the opposite behavior, with negatively 
skewed distributions in the inner Galaxy and positively 
skewed distributions in the outer disk. For this paper, 
we did not correct for the APOGEE selection function. 
In the future, we plan to do a more detailed comparison 
between APOGEE observations and the simulation from 
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Minchev et al. (2013), in which the selection function is 
taken into account. 


4.2. Radial Mixing 

Simple chemical evolution models (closed or leaky box) 
are unable to produce the positively skewed MDFs that 
we observe in t he outer disk. Mo dels th at include ra¬ 
dial mixing (e.g..Schonrich & Binney|2009| are able to at 

an-shaped ME 


least produce a more Gaussian-shaped 1V1DF across the 
disk. The fraction of stars that undergo radial migra¬ 
tion is difficult to predict from first principles because it 
depends in detail on spiral structure, bar perturbations, 
and perturbations b y and merge r s with sate lli tes (e.g., 


Roskar et al. (2008); Bird et al. (2012, 2013); Minchev 
et al.| (|2013|)). 

To test the effects of blurring and churning on our 
observed MDFs, we create a simple model of the MDF 
across the disk. We assume that the intrinsic shape of 
the MDF is uniform across the disk, and that the ob¬ 
served change in skewness with radius is due to mixing 
of populations from different initial birth radii. The disk 
is mo deled as a Dehnen distribution function (Dehnen 
1999) with a velocity dispersion of 31.4 km s -i , a ra¬ 
dial scale length of 3 kpc, and a flat radial velocity dis¬ 
persion profile. These distribution function parameters 
adequ ately fit the kinematics of the main APOGEE sam¬ 
ple (Bovy et al. 2012d). We model the initial MDF as 
a skew-normal "distribution with a peak at +0.4 dex in 
the inner Galaxy, a dispersion of 0.1 dex, and a skewness 
of —4; we assume a radial gradient of —0.1 dex kpc -1 
to shift the peak of the MDFs as a function of radius, 
keeping the dispersion and skewness fixed. We then de¬ 
termine the distribution of guiding radii with blurring or 
churning to determine the effect on the observed MDFs. 


4.2.1. Blurring 



[Fe/H] 


Figure 8. The MDF as a function of R for our simple blurring 
model. The dashed lines show the initial MDFs (which the stars 
on circular orbits would have) and the solid lines represent the 
MDFs including the effect of blurring. The skewness of the blurred 
MDFs are indicated. The magnitude of the skewness diminishes 
with radius, but it does not change sign. 


To determine the effect of blurring on the observed 
MDF, we use the model as described above and com¬ 
pute the distribution of guiding radii R g due to blurring. 
Assuming a flat rotation curve and axisymmetry, R g V c 


is equal to RV r Q t, where Y ro t is the rotational velocity 
in the Galactocentric frame. The blurring distribution 
Pb(R g \R) is given by 

p b (R g \R) ~ p(V Iot = ^V C \R). 

The probability p(V rot \R) can be evaluated using the as¬ 
sumed Dehnen distribution function. The resulting MDF 
is 


P([Fe/H]|i?) = J dR g p([Fe/H\\R g )pb(R g \R) 

In Figure |8] we compare the initial MDF and the MDF 
with the effects of blurring included. While blurring does 
reduce the observed skewness of the MDFs, the MDFs 
are still negatively skewed at all radii. This model is 
simplistic and it is possible that our underlying assump¬ 
tion regarding the intrinsic shape of the MDF may not 
be correct. However, because the intrinsic MDF is un¬ 
likely to have positive skewness anywhere in the Galaxy, 
it appears that blurring alone is unable to reproduce 
the change in sign of th e MDF sk ew ness se en in the 
APOGEE observations. |Snaith et al.| (|2014|) find that 
solar neighborhood observations can be adequately ex¬ 
plained by the blurring of inner and outer disk popula¬ 
tions, but we conclude that such a model cannot explain 
the full trends with Galactocentric radius measured by 
APOGEE. 


4.2.2. Radial Migration 

We expand our simple model to include churning to 
determine if radial migration is better able to reproduce 
the observations. We model the effect of churning on the 
guiding radii by assuming a diffusion of initial guiding 
radii R g ^ to final guiding radii R g j , for stars of age (r), 
given by 


P(Rgj\ R g,i’ T ) = 

Af {gigj\Rg^i, 0.01 + 0.2 T R, 


0 -(Rg,i- 8 kpc) 2 /16 kpc' 


')• 


where J\f(-\m,V) is a Gaussian with mean m and vari¬ 
ance V. The spread increases as the square root of time 
and the largest spread in guiding radius occurs around 
8 kpc. Analysis of numerical simulations has found that 
migration ca n be accurately described as a Gaussian dif¬ 


fusion (e.g., Brunetti et al. 2011 Kubryk et al. 2014 
Yera-Ciro et al. 2014p ; we use the simple analytic form 


above"to+tpproximate the effect of churning seen in these 
more realistic simulations. 

To obtain the churning distribution p c (R g \R, r) of ini¬ 
tial R g at a given radius R and age r, we need to convolve 
p(R g j\R g ^i : r) with the blurring distribution pb(R g \R) 
above 

Pc(Rg\R,r) = J dRgjPb(Rg,f l-R) P{Rg,f \Rgi T ) 5 (1) 


where we have used the fact that churning leaves the total 
surface density approximately unchanged and therefore 
that p(R g \Rg if ,r) ~ p(R g , f \R g , t). 

In order to determine the MDF due to churning, we 
must also integrate over the age distribution at a given 
radius and therefore need (a) the metallicity as a function 
of age at a given initial radius and (b) the age distribution 
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at every R. For the former, we assume that metallicity 
increases logarithmically with age as a function of ra¬ 
dius, starting at [Fe/H]= —0.9 and up to the peak of the 
skew-normal intrinsic MDF plus its dispersion; this rela¬ 
tion is shown for a few radii in Figure [9] We approximate 
the final age distribution at each radius as the initial age 
distribution, which is approximately the case in a more 
detailed calculation that takes the effects of churning into 
account. The initial age distribution p(r\R) at each ra¬ 
dius is simply a consequence of the assumed initial MDF 
and initial age-metallicity relation. 

We obtain the churning MDF as 

p([Fe/H]|ii) = J dRgp c (R g \R,r) p(r\R)/ |d[Fe/H]/dr| , 

where the age r is a function of R g and [Fe/H]. We eval¬ 
uate this expression at different radii and obtain MDFs 
that include the effects of both blurring and churning, 
displayed in Figure [lOj With the addition of churning, 
we are able to reproduce the change in skewness observed 
in the MDFs across the plane of the disk and in partic¬ 
ular the change in sign around R = 9 kpc. Stars need to 
migrate at least 6 kpc to the outer disk and at least 3 kpc 
around the solar neighborhood to produce the observed 
change in skewness. 

The model used in this section is a very simple model, 
and the predictions from the model do not match our 
observations perfectly; more detailed and realistic mod¬ 
eling is required that consistently takes into account the 
chemo-dynamical evolution of the disk. However, these 
tests demonst rate that blurring alone, as suggested by 
laith et al. (2014|), is unable to reproduce our obser¬ 
vations, and that the addition of churning to our model 
yields significantly better agreement with the observed 
MDF, in particular the change in skewness with radius. 
Therefore, the APOGEE MDFs indicate that migration 
is of global importance in the evolution of the disk. 

In principle, the gas can undergo radial mixing as well, 
causing enrichment (or dilution) of the ISM at other lo¬ 
cations in the Galaxy through processes such as outflows 
and Galactic fountains. However, in this case the gas 
must migrate before it forms stars, so the timescales are 
much shorter. Additionally, many processes that would 
cause gas to migrate, such as non-axisymmetric pertur¬ 
bations, will also induce stellar migration and can there¬ 
fore not be decoupled from the stars. It is likely that a 
combination of both gas and stellar migration is required 
to reproduce our observations in chemo-dynamical mod¬ 
els for the Milky Way. 

Finally, we can also calculate the predicted age- 
metallicity relation at different locations throughout the 
disk from this model: 

p(r,[Fe/U\\R) = p(Rg\R,T)p(r\R)/ \d[Fe/H\/dR g \ , 

where R g is the initial guiding radius corresponding to r 
and [Fe/H]. The age-metallicity relation from the model 
is shown in Figure [9] This relation is in qualitative agree¬ 
ment with observations of the age-m etallicity r elation ¬ 


ship of the solar neighborhood (e.g., Nordstrom et al. 


2004 Haywood et al. 


2013 


Bergemann et al. 


2014jrf 
f metallic- 


4.3. Vertical Structure: Heating vs. Cooling 

Our measurements confirm the well established trend 
of enhanced [a/Fe] in stars at greater distances from the 
plane, and they show that this trend with \z\ holds over 
radii 3 < R < 11 kpc. At still larger radii, the fraction 
of high-[a/Fe] stars is small at all values of |z|. The 
fact that stars with thin disk chemistry can be found at 
large distances about the plane in the outer disk can be 
explained by viewing the disk as composed of embedded 
flared mono -age populations, as suggested by Minchev 
eFaTl([M5). 

The trend of [a/Fe] with \z\ is particularly striking 
in the 3 — 5 kpc annulus, where the stars lie along the 
sequence expected for a single evolutionary track but the 
dominant locus of stars shifts from low-[Fe/H] , high- 
[a/Fe] at 1 < \z\ < 2 kpc to high-[Fe/H] , low-[a/Fe] at 
0 < M < 0.5 kpc. Along an evolutionary track, both 
[Fe/H] and [a/Fe] are proxies for age. 

Heating of stellar populations by encounters with 
molecular clouds, spiral arms, or other perturbations will 
naturally increase the fraction of older stars at greater 
heights above the plane simply because they have more 
time to experience heating. However, even more than 
previous studies of the solar neighborhood, the transition 
that we find is remarkably sharp, with almost no low- 
[a/Fe] stars at |z| > 1 kpc and a completely dominant 
population of low-[a/Fe] stars in the midplane. Explain¬ 
ing this sharp dichotomy is a challenge for any model that 
relies on continuous heating of an initially cold stellar 
population. One alternative is a discrete heating event 
associated with a satellite merger or other dynamical en¬ 
counter that occurred before SNIa enrichment produced 
decreasing [a/Fe] ratios. Another alternative is “upside 
down” evolution, in which the scale height of the star¬ 
forming gas layer decreases with time, in combination 
with flaring of mono-age populations with radius (e.g., 
Minchev et al.| (|2015|) ), as the level of turbulence associ¬ 
ated with vigorous star formation decreases (e.g., |Bour 


our model, the outer disk has a wide range o 
ities at any age due to radial m igration, in agreement 
with models from Minchev et al. (2014). 


naud et al.|2009| |Krumholz et al.|2012||Bird et al.|2U13|). 

In this scenario, the absence of low-[a/Fe] stars far from 

the plane implies that the timescale for the vertical com¬ 
pression of the disk must be comparable to the 1 — 2 Gyr 
timescale of SNIa enrichment. This timescale appears at 
least roughly consistent with the pre dictions of cosmo - 
logical simulations (e.g., Figure 18 of Bird et al. 2013). 
More generally, the measurements presented here of the 
spatial dependence of [a/Fe] -[Fe/H] distributions and 
MDFs provide a multitude of stringent tests for models 
of the formation and the radial and vertical evolution of 
the Milky Way disk. 

5. CONCLUSIONS 

The solar neighborhood MDF has proven itself a 
linchpin of Galactic astronomy, enabling major ad¬ 
vances in our understa nding of chemical evolution (e.g., 
van den Bergh 1962 ) in conjunct ion with stellar dy¬ 
namics (|Sch6nrich & Binney 2009). Galaxy formation 
models must ultimately reproduce the empirical MDF 
as well (Larson 1998). In this paper, we report the 
MDF and alpha-abundance as measured by 69,919 red 
giant stars observed by APOGEE across much of the 
disk (3 < R < 15 kpc, 0 < |z| < 2 kpc). Our simple 
dynamical model reveals the exciting prospect that the 
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Figure 9. The age-metallicity relation at different radii for our simple model. The blue line shows the metallicity of the ISM as a function 
of age at different radii, while the solid black lines denote 1,2, and 3 cr of the distribution of stars in the age-metallicity plane respectively; 
the solid red line gives the mean metallicity as a function of age. The solar neighborhood and the outer disk both display a wide range of 
metallicities for all but the oldest stars. 



[Fe/H] 


Figure 10. The MDF as a function of R with the inclusion of 
blurring and churning. The dashed lines show the initial MDF and 
the solid lines display the MDF including the effects of churning 
and blurring. The redistribution of guiding radii due to churning 
is able to significantly change the skewness compared to the ini¬ 
tial MDFs and can explain the changes in skewness and its sign 
observed in the APOGEE sample. 

detailed shape of the MDF is likely a function of the dy¬ 
namical history of the Galaxy. Our conclusions are as 
follows: 

• The inner and outer disk have very different stellar 
distributions in the [a/Fe] vs. [Fe/H] plane. The 
inner disk is well characterized by a single sequence, 
starting at low metallicities and high-[a/Fe] , and 
ending at solar- [a/Fe] abundances with [Fe/H] ^ 
+0.5, while the outer disk lacks high-[<a/Fe] stars 
and is comprised of primarily solar-[a/Fe] stars. 

• The scale-height of the inner Galaxy decreased with 
time, as the (older) metal-poor high-[a/Fe] stars 
have large vertical scale-heights, while (younger) 
metal-rich solar-[a/Fe] populations are confined to 
the midplane. 

• The peak metallicity and skewness of the MDF is 
a function of location within the Galaxy: close to 
the plane, the inner disk has a super-solar metallic¬ 
ity peak with a negative skewness, while the outer 


disk is peaked at sub-solar metallicities and has a 
positive skewness. 

• Models of the MDF as a function of R that in¬ 
clude blurring are unable to reproduce the observed 
change in skewness of the MDF with location. 

• Models with migration included match our obser¬ 
vations of the change in skewness; migration is 
likely to be an important mechanism in the ob¬ 
served structure of the disk. 
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APPENDIX 

A. POTENTIAL METALLICITY BIASES 

We consider several potential sources of bias in our observed metallicities and their distributions: 

• Bias from sample selection: color cut and magnitude limits 

• How well a selection of giants is able to reproduce the underlying MDF of the sampled population 

• ASPCAP spectral grid edges: spectral library has no results for T e ff < 3500 K and log# < 0 

• Distance Errors 


A.l. Correcting For Biases From Sample and ASPCAP Selection 

We use TRILEGAL simulations (Girardi et al.| 2005) of stellar counts across the Galaxy to determine the magnitude 
of any sample selection effects on our results. These simulations provide an opportunity to test how well our targeted 
sample is able to trace the underlying population of giants as a function of location within the disk. Note that giants 
are not a perfect representation of the underlying populations of stars, having a bias against older (generally more 
metal-poor) stars, with the relative population sampling being roughly proportional to r -0 6 . However, correcting for 
this bias is difficult, and depends strongly on the star formation histories at different locations throughout the disk. 
For this work, we do not correct for the uncertain, non-uniform sampling of giants and therefore our MDFs are slightly 
biased against older populations. 

The simulations were run through our targeting software, to obtain an APOGEE-like sample through all of the 
lines of sight observed in the DR12 dataset. These simulations do not have a radial metallicity gradient, so they are 
not a perfect representation of our observations, but are still useful for our testing purposes. Because of the lack 
of a metallicity gradient, the inner disk (3 < R < 7 kpc) is significantly more metal-poor in the simulations than is 
observed, which means the stars making up the stellar populations in the simulation have higher effective temperatures 
than would be observed in these locations. Because of this effect, we also include the bulge regions from the simulations 
(0 < R < 3 kpc) as it is more metal-rich, to aid in our testing of the effects of the ASPCAP temperature grid edge on 
the obser ved metallicity distributions. The underlying MDF is compared to the observed MDF as a function of radius 
in Figure [Alj The observed MDF is similar to the underlying MDF based on our sample sele ctio n - we do not have 
large systematics due to our field selection, magnitude limits, or color cut (panel (A), Figure Al). We also compare 
the ages of the stars we observe to those of the underlying giant populations; the age distribution for our targ eted 
sample of giants is in general older than the underlying sample, but the differences are not dramatic (Figure A2). 

However, in the inner Galaxy (0 < R < 3 kpc) there is a significa nt b ias against metal-rich stars due to the 
temperature edge in the ASPCAP spectral libraries (panel (B), Figure Al). This result is expected, as the more 
distant and more metal-rich stars are cooler. The magnitude of this bias on our observed APOGEE sample depends 
on the number of metal-rich stars present at larger distances, but it is likely to be significant for stars with [Fe/H] 
>0.1 based on these simulations. There are several thousand stars in the DR12 catalog that are at the temperature 
grid edge; most are in fields towards the inner disk and Galactic center, and are likely the most metal-rich stars in the 
sample. It is therefore likely that there is a significant bias against metal-rich stars for these fields. 

To mitigate this bias, we apply a cut in surface gravity on our sample. Marching down the giant branch to higher 
surface gravities lowers our sample size, especially at larger distances, but reduces the impact of the temperature grid 
edge on our sample (demonstrated in Figure [TJ. The surface gravity cut removes the metal-poor, luminous giants 
that are warmer than the grid edge than the similarly lower surface gravity metal-rich stars, which are removed from 
the sample due to the temperature restrictions. A cut of lo g q >1.0 was most effective at mitigating the bias due to 
the temperature grid edge (panels (C) and (D), Figure [AT] ). There is still a slight bias against metal-rich stars even 
with this surface gravity cut, but it is significantly less than if no restriction on surface gravity was imposed. It is 
possible that the observed MDFs for the very inner disk (3 < R < 5 kpc) are slightly more metal-poor than in reality, 
but the effect is likely not large; the difference in mean metallicity of the MDF from the input simulations and our 
observed simulations is ^0.1 dex after the surface gravity restriction, compared to > 0.3 dex due to the effects of the 
temperature cutoff. 
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Figure Al. The observed mock MDF of giants (red) compared to the intrinsic MDF of giants (black) as a function of R. (A): Only 
targeting criteria (field placement, color cuts, magnitude limits) applied to the input simulations. (B): The observed MDF with the 
temperature restrictions imposed by the ASPCAP grid edge. (C): The observed MDF with a restricted surface gravity range. (D): The 
observed MDF with both the temperature and surface gravity restrictions. There is little difference between panels (C) and (D), the 
effects of the temperature grid-edge on the observed MDF has largely been mitigated due to the restriction in surface gravity, leading to a 
significantly less biased MDF. 
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Figure A2. The intrinsic age distribution of giants in the TRILEG AL simulation (black) compared to the mock observed age distribution 
(red). The sample selection criteria does not induce a large bias in the observed sample, although the targeted sample is slightly older in 
general than the underlying age distribution of giants in the simulation. 


A.2. Distance Errors 

Distances to individual stars are accurate to about 15-20%. With nearly 70,000 stars in the sample, there are 
non-negligible amounts of 2-3cr errors in the distances, which could potentially influence the observed metallicity 
distributions. In particular, the largest distance bins (inner and outer disk) will suffer more from errors in distance 
calculations, as 20% errors in distance at these locations is roughly the radial extent of the bin. We use TRILEG AL 
simulations to test the effect of distance errors on our sample. 

A radial gradient of -0.075 dex kpc -1 is added to the default MDF of the simulation, such that the inner disk is 
more metal-rich than the outer disk. Distance errors are modeled using a normal distribution with a dispersion of 
20% and mean of 0. The MDF of the targete d sim ulations (the “true” mock observed MDF) is compared to the mock 
MDF with distance errors included in Figure |A3| The errors in distance have a negligible effect on the MDF, except 
for potentially in the very outer disk (R > 13 kpc). This area of the Galaxy has the fewest stars per R,\z\ bins as 
well, so any scatter from other regions of the disk will have a larger impact on the observed MDF. The change in the 
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[Fe/H] 


Figure A3. The mock targeted MDF (black) versus the mock observed MDF after applying distance errors to the TRILEGAL simulation 
(red). The difference between the two is negligible except for the very outer disk (R > 13 kpc). 
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Figure Bl. The distribution of stars in the [a/Fe] vs. [Fe/H] plane color coded by T e ff for stars with \z\ < 0.5 kpc. There is evidence for 
temperature trends in the [a/Fe] abundance in the inner Galaxy (3 < R < 5 kpc). This figure is down-sampled and shows only 15% of the 
datapoints for R > 5 kpc. 
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Figure B2. Similar to above, but with the surface gravity restrictions removed. The [a/Fe] abundance trends with temperature are clear 
with the inclusion of the lower surface gravity stars. Cooler stars have lower [a/Fe] abundances by ~ 0.1 dex. 


MDF due to distance uncertainties even at these larger radii is still small, however, less than ~ 0.1 dex, so the overall 
effect on our results as a whole due to distance errors is negligible. 

B. SYSTEMATICS IN [a/Fe] WITH TEMPERATURE 

Th ere are temperature trends in our [a/Fe] abundance measurements. Calibrations performed by Holtzman et al. 
(2015) remove most of these trends, but they lack calibrators for cool stars (T e ff < 4000 K) in most clusters. For 
stars cooler than this value, there are still temperature trends in [a/Fe] . The trends in [a/Fe] with temperature likely 
reflect that different elements contribute to the overall [a/Fe] determination at different temperatures. 

The surface gravity restrictions employed to mitigate the metallicity bias due to the temperature cutoff in the spectral 
libr aries removes most stars that are affected by errors in the [a/Fe] determination, but some issues remain (Figure 
Bl). The warmer stars have higher [a/Fe] abundances, while the cooler stars have slightly lower [a/Fe] abundances. 
'The overall spread due to temperature is small, ~ 0.1 dex, but is noticeable. It is also possible that the bowl shape 
of the low-[a/Fe] sequence observed in the solar neighborhood is due to similar temperature effects: the cooler stars 
again have lower-[a/Fe] abundances, blurring the overall shape of the low-[a/Fe] sequence. 

If the restrictions in su rfac e gravity are not enforced, the appearance of the [a/Fe] vs. [Fe/H] plane is different in 
the inner Galaxy (Figure |B2| ), and the temperature effects on [a/Fe] abundance become more apparent. These stars 
were removed due to the restriction in surface gravity, but are also the coolest stars in the sample that have the most 
dramatic errors in the [a/Fe] determination. The fact that there is not a higher surface gravity component to this 
population is a clue that these stars may have a poor [a/Fe] determination. We also meas ure the [a/Fe] abundance 
by averaging over the abundances of the individual a elements (O, Mg, Si, S, Ca, Figure |B3|) , and find that these 
stars show a similar temperature gradient compared to the overall [a/Fe] measurement, with the cooler stars having 
lower-[a/Fe] abundances by ~ 0.1 dex at the same [Fe/H] . It is likely that there are errors in the [a/Fe] abundances for 
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Figure B3. Simila r to above, but with the [a/Fe] measured by averaging the abundances for the individual a elements. We find similar 
trends as in Figure. |B2| above, with the cooler stars having lower-[a/Fe] abundances than the warmer stars at the same metallicities. 

these cool low-surface gravity stars, that are removed by our corrections for the metallicity bias due to the ASP CAP 
temperature grid edge. It is possible, however, that these abundances are accurate and that there is a second sequence 
of stars present in the inner disk, this will affect the interpretation of our results for the inner Galaxy. The [a/Fe] 
abundances for these cooler stars are higher than that of the solar- [a/Fe] populations observed in the rest of the disk, 
but slightly lower than the high-[<a/Fe] sequence, and if real are a potentially unique stellar population. 











